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ABSTRACT 

The mathematical model of process of extraction of 
an oilseed material in the units with critical carbonic 
gas was developed. On the basis of this model 
researches are conducted and schedules of change of 
concentration of oil in firm both liquid phases and 
decrease of concentration of oil in quasilayers are 
gained at different coefficients of diffusion, and also 
resulted comparison of the gained results to 
experimental data. 

The main property of technological processes is that 
the components of the complex's events have a 
deterministic-stochastic nature, this nature is 
manifested in mass-exchange and heat-mass exchange 
processes and in the restoration of stochastic 
properties of the hydrodynamic medium in chemical 
transformations. The components of these phases are 
characterized by a random interaction of the 
constituent components or by the boundary conditions 
of the geometrical characteristics of the apparatus 
(productive orientation of the inter phase boundary of 
the moving media). Similarly, different systems are 
characterized by a complex interaction of phases and 
component and parts of components, as a result of 
which it is impossible to study these systems from the 
standpoint of the laws of conservation or transfer of a 
classical deterministic substance. 
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To solve a similar problem, it is necessary to apply 
methods of system thinking, object analysis, 
mathematical modeling and finding the optimal 
solution. The basis of this strategy is the adoption of 
quasi-objects as a complex interacting hierarchical 
system, the conduct of qualitative analysis of 
structures, evaluation and definition of parameters, 
their mutual influence, in particular, in the form of a 
mathematical description of the process [1,2,3]. 

The process of extraction of substances with 
supercritical carbon dioxide (SC CO 2 ) has principal 
advantages. The properties of CO 2 in the supercritical 
state are intermediate between its properties in the gas 
and liquid phases. So, SC CO 2 , like all supercritical 
gases, has a high density, close to the liquid, and low 
viscosity as gases. The density of SC CO 2 differs 
slightly from the density of the liquid, i.e. SC CO 2 has 
a significant capacity, the viscosity is close to the 
viscosity of the gas, and the diffusion coefficient has 
an intermediate value, which significantly accelerates 
mass transfer processes. This combination of 
properties makes it extremely effective to use SC CO 2 
as substitutes for organic solvents in laboratory and 
industrial processes [4,5,6], 

The process of extraction of plant matter can be 
divided into 3 stages [7]: 1) "internal diffusion", 
which includes all the phenomena of transport of 


@ IJTSRD | Available Online @ www.ijtsrd.com | Volume - 1 | Issue - 6 | Sep-Oct 2017 


Page: 50 







International Journal of Trend in Scientific Research and Development (IJTSRD) ISSN: 2456-6470 


substances within the raw material particles (solvent 
penetration into the pores of the particles of plant raw 
materials, dissolution of the target component (s), 
transfer of the extractable matter within the particle of 
the plant raw material to the phase interface, 2) the 
transfer of matter within the boundaries of the directly 
diffusion boundary layer; 3) the transfer of the 
extracted substance, the moving extractant from the 
interface and the distribution of it over the entire mass 
of the extractant (convective diffusion). 

In mathematical modeling of the extraction process, 
the abstract representation of the process is of decisive 
importance. Below, a mathematical modellingof mass 
transfer in the quasi-layer in the “solid particle-gas- 
liquid” system is presented. 

To do this, you need to enter some restrictions and 
assumptions: 

> particles of oil-containing material (cake) have an 
ideal ball shape; 

> the particle of the cake is divided into quasi-layer 
and these layers have the same thickness; 


> in each layer of the particle the oil content is 
distributed evenly; 

> oil is only in the pores of the particle; the oil 
content of the material is equal to its porosity; 

> For each quasi-layer, the previous layer is 
"external", the next layer is internal, respectively. 

Keywords: system; analysis; modeling; quasi-layer; 
diffusion; extraction; extraction unit; oilseed 
material; carbon dioxide, supercritical condition 

Modeling of the extraction process in the quasi¬ 
layer of the material particle. Taking into account 
above assumptions, we will form a mathematical 
model of the mass change in the internal and average 
quasi-particles of the particle, as well as in quasi-layer 
directly in contact with the solvent. As is known, in 
matter quasi- media the substance passes to the 
surface of contact with the solvent through molecular 
diffusion. For example, fig. 1 shows the input and 
output parameters of the extraction process in the 
medium quasi-particles of the cake particle 
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Fig.l. Input and output parameters of the extraction process in medium quasi-layers particles of a cake 

In fig. 1, the quasi-layer is numbered consequently starting from the outer layers. 

The equation of material balance in average is a quasi-layer material particle: 


oU3 _ Q _Q 

com 3 w con 3 


( 1 ) 


Where: G CO m3 - oil in the quasisof the particle cake, (m 3 /s); G co „3 - oil consumption through a quasis particle 
cake, (m /s). 


As is known, the amount of oil M 0 n 3 passing through the medium quasi-layer the outer layer is the product 
number micelleM,„ t j, this last layer of the concentration of oil in the micelles,,//}: 


M = M , • a 

oil 3 me 3 oil 3 


( 2 ) 


Then: 


djMjnc 3-0 


dt 


— G — G 

W com 3 W con 3 


(3) 


Next, we analyze the mathematical expressions describing the changes in the concentration of oil in the micelle: 


da 3 

dt 


m 3 


(G 


com 3 


^con 3 ) 


(4) 


Mass micelle/n.?, located quasi-layer average particle equal to the product volume F^and density p mC 3 this 
micelle: 


m 3= V ,nc3-P f 


mc 3 


(5) 


The volume of the micelle V mC 3 located quasi-layer average particle equal to the product of the difference 
between the volumes of the outer layer V out3 a.nd the inner layer V inn4 and porosity (oil) cakes z: 

V m c3 = { V ou,3 - V inn 4 )- Z ( 6 ) 

Since the layers have the shape of a sphere, their volumes aredetermined by the following formulas: 


V t =--n-r \ 3 

out 3 ^ out 3 


V =-.^. r 3 

inn 4 ^ inn 4 


(7) 


where: r ou t 3 and r„,„^-radii of the outer and inner quasi-layer of the cake, respectively. 
We determine the radii of the outer and inner quasi-layer of the cake particle: 
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Micelle density: 


=-•(«-* + !) 

n 

r =-•(«-/) 
n 


( 8 ) 


P mC 3=P •0- a i ) + Po,r a o i « ( 9 ) 


Where p and p 0 u are the density of the solvent and oil, respectively, (kg/m ). 

Oil flow passing from the average quasi-layerthe outer layer is the product of the mass transfer coefficient, and 
the mass transfer area concentration difference (driving force of the process): 

G con3= k -Fo Ut -{a i n tB -a t2 ) ( 10 ) 

G cona= k - F u\ a , 4-0 (H) 


Mass transfer coefficient is directly proportional to the diffusion coefficient and the density and inversely 
proportional to the micelle layer thickness, i.e.: 


k = 



( 12 ) 


The thickness of the quasi-particle of the cake particle <5is determined by the ratio of the particle radius R to the 
number of layers n: 



n 


(13) 


The working areas of the outer and inner quasi-layer of the cake particle are determined: 


F 


out 3 


F 


inn 4 


A-71 

A-71 



(14) 
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Fen Vi bail 


Fig.2. Computer model of an average quasi-layer of a particle of an oilseed material 

Figure 2 shows a computer model of the mean quasi-layer of a cake particle, which consists of a block of the 
radius of the inner layer Fcnri, and its dependent blocks: the outer area of the inner layer Fen Fi, the volume 
block of the inner layer Fen Vi bail and the block of oil from the inner layer FenG com . 

A mathematical description of the average concentration of oil in a micelle, located in the quasis cake 
carried out. In order to increase the adequacy of the results of the mathematical model and the experimental 
results in the calculations, we assume that the average concentration of the oil of the oil cake leaving the quasi¬ 
layer is equal to the concentration of the oil leaving the particle to the solvent, i.e. ao = a ave . 

The average concentration of the oil of the oil cake leaving the quasisis determined by the following 
mathematical expression: 


5 5 

a n — a = V m a /'Y m . (15) 

0 ave i i i V / 

i=l i =1 

i.e. the average concentration of the oil of the oil cake leaving the quasi-layer is equal to the ratio of the 
products of mass and oil concentration in each quasis to the sum of the oil mass in these quasis: 

ao=a ave =(miai+ m 2 d 2 + msa 3 + ni 4 a 4 +msa 5 )/ mi+ ni 2 + m 3 + ni 4 + m 3 (16) 

where: mi, m 2 , m 3 , m 4 , m 3 - the oil mass in each quasi-layer, respectively; aj, a 2 , 013 , CI 4 , as - the oil concentration 
in each quasi-layer, respectively. 
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Fig.3.The average oil concentration in the particle quasi-layer during the extraction of the oilseed 

material by supercritical CO 2 


Figure 3 shows a mathematical model for calculating 
the average concentration by creating a common 
block of oil concentration emerging from the quasi¬ 
layer, the concentration of oil entering the next layer, 
and the mass of the oil in the quasi-layer. The average 
oil concentration in quasi-shells in a computer is 

expressed in the form FcnG com -u (1) . u (10). For 

the given model, in the form of input values, oil 
concentrations are taken in the form of an oval 
emerging from each quasis and a mass corresponding 
to these layers. The value of the average oil 


concentration in the cake particle leaving the Fen M 
block in the form of an oval is adopted as an input 
value for the next computer model. 

Modelling of the process of extraction of cake in a 
unit with critical carbon dioxide. As is known, on 
the basis of the material balance of extraction, it can 
be said that during extraction, the oil passing from the 
composition of the material to the solvent changes its 
consumption and the consumption of the solid does 
not change. Taking this into account we modelling the 
process in one level of the unit. 


G 


com * 1 



Fig.4. The scheme of the extraction process in the suspended state of the cake 
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Figure 4 shows incoming and outgoing indicators to the apparatus as an extraction object, the following: 
incoming cake feed, initial oil content of cake, cake consumption, oil content of the cake, incoming solvent 
consumption and concentration, and the concentration of the micelle. 

The materialbalance (fig.5,6),which confirms the equality of the difference between the consumption of the 
solid phase (oilcake cake) and the oil consumption in its composition and the difference between the 
consumption of cake and the residual quantity of oil in the cake, expresses the invariance of the consumption of 
the solid phase during the extraction of the cake of the oilseed material. 

This material balance is as follows: 


Go- Go'ao -Gc-Gc'Ci (17) 

where: Go - cake consumption in the apparatus; ao- is the initial oil content of the cake, i.e. the average 
concentration of oil in the micelle ( a ave ) leaving the layers of the cake particle; G c - is the consumption of the 
cake from the apparatus; a-oil content of the cake leaving the machine. 

In (17), we derive from the brackets the same variables and obtain: 

Go-(l-ao) = G c (l-a) (18) 

Consumption of cake from the machine: 

G, = G 0 (19) 

1 - a 

The consumption of oil, which has passed into the liquid phase, i.e. in the solvent isdetermined by the 
difference in oil consumption when entering the machine and the oil consumption when the machine leaves: 

G 0 u=Go'Cio-G-a (20) 

or, if G expresses through Go, then we get: 

G„„ (21) 

1 - a 

In general, the oil consumption from oilseed oil cake in the solvent, which depends on GO, ao, a, is determined 
by the following formula: 


G OII = G 0 • (a 0 -1—— • a) (21) 

1 - a 

According to the material balance in the "solid-liquid" phase system in the first level of the apparatus, we 
determine the change in the concentration of oil in the solvent in time: 

^ = ^-(Ly„-Ly + GJ ( 22 ) 
dr Vp 

where: L, V and p- are the flow rate, volume and density of the liquid phase, respectively; >vr is the initial 
concentration of oil in the solvent; y- is the output concentration of the oil in the micelle; G 0 n- is the 
consumption of oil that has passed from the cake to the solvent. 
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From the equality of solvent consumption at the entrance and exit from the first level of the apparatus, we 
obtain: 


Lo-(l-yo) =L(l-y) (23) 

where: Z 0 -consumption of solvent at the entrance to the apparatus; Z-consumption of the micelle from the first 
level of the apparatus. 

Hence the expenditure of the missile from the first level of the apparatus is equal to: 

L = L 0 (24) 

i -y 


According to equations (22) and (24), we obtain: 


dy 

—' ™ = V To 
dr 


L»- —y + G oll 


i -y 


or 


±. m= L o(yo -^ y) + G 

dr 1 — y 


oil 


(25) 


(26) 


Carrying out some arithmetic operations in expression (26), we determine the concentration of the micelle 
when leaving the level y: 


T = 


Vo + G, 


oil 


G oa + L 


(27) 



Fig.5.The computer view of the process model blocks in the working chamber of the extraction apparatus 
of the cake of oilseed material with supercritical carbon dioxide 
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Fig.6.The computer view of the process model block in the solid phase in the working chamber of the oil 

cake extraction unit with supercritical carbon dioxide 


In the computer model (fig. 7), the oil consumption from oilseed oil cake to the solvent can be represented 
graphically in the Fen M block using scope a2. This block enters as an input parameter into the Fcndy / dt block 
and in this block of the computer model the concentration of the micelle in the first level of the apparatus and 
the residual oil content of the cake are reflected. 



Fig.7.Computer view of the model block of the process in the liquid-gas phase in the working chamber of 

the oil cake extraction unit with supercritical carbon dioxide 



Siklon 3 


Fig.8. Expression of the process model in the working chamber 


Figure 8 shows the general model of the extraction 
process of the oil cake in the extraction chamber. To 
work with this model, we have formed a model for 
introducing initial values of the extraction parameters, 


i.e. critical pressure of carbon dioxide as an 
extractant, its temperature, etc. Several experiments 
were carried out with the model developed for 
different diffusion coefficients. After the introduction 
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of the initial parameters, when pressing the 
"Calculation" button, it is possible to obtain a graph of 
the change in the oil concentration in the solid and 
liquid phases of the material (fig. 9a) and the decrease 
in the oil concentration in the quasilocks of the cake 
particle having certain diameter values (fig. 9b). As 
can be seen from fig.(fig. 9b), the concentration of oil 
in the outer layer decreases faster, i.e. here oil is 
relatively quickly removed to the outer layer and oil is 
supplied from the lower layer. The intensity of mass 
transfer depends on the value of the diffusion 
coefficient. In our case, the value of the diffusion 
coefficient is chosen to be Z)=4.5T0' 12 m 2 /s, since in 
the case of pre-crushed and processed oil-based raw 
materials, such an index of the diffusion coefficient 
corresponds, according to [8], The extraction time is 
600 seconds. At the same time, oilseed oil content 
decreases from 15% to 3%. 

When the raw material is not processed, as shown in 
the sources [8], the diffusion coefficient is an order of 


es* 5 11 



a 


magnitude lower, D=4.5T0' 13 m 2 /s. Figure 10 shows 
the graphs of the change in oil in the solid and liquid 
phases (a) and the decrease in the concentration of oil 
in the quasi-layer (b) with the above-mentioned 
diffusion coefficient obtained in the mathematical 
model. The nature of the process is similar, but in this 
case, to reduce the oil content of oilseed oil cake from 
15% to 3%, 4000 s time is required, as the process of 
diffusion of substances from undisturbed raw material 
cells is slow, so the extraction rate in this case is slow. 
The diffusion coefficient, which depends on the 
internal structure of the plant cell, changes in the 
content of substances in time, changes in the cellular 
structure and concentration of diffusing substances 
during extraction, on the properties of the extragent 
and the extractable substance, and also on the 
temperature and hydrodynamic conditions of the 
process. 


e* 5 1 1 



6 


Fig.9. Change in the concentration of oil in the solid and liquid phases with a diffusion coefficient D = 
4.5*10' 12 m 2 / s (a) and a decrease in the oil concentration in the quasi-layer (b). 
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Fig.10. Change in the concentration of oil in the solid and liquid phases with a diffusion coefficient 
D = 4.5*10' 13 m 2 /s (a) and a decrease in the oil concentration in the quasi-layer (b). 


Thus, using of the developed mathematical model, a 
theoretical picture of the process was considered. 

In order to compare the theoretical experiments with 
the experimental ones and to determine the rate of 
extraction of the lipid complex from the seeds of the 
grapes, a series of experiments was carried out on the 
following parameters: grinding of the grape seeds to a 
particle size of 1-2 mm and lapping on the rolls. The 
thickness of the petal was 0.2 mm. Extraction mode: 
temperature 20°C, pressure - 6,5MPa, extraction time 
- 600 s. The content of lipids in the feedstock was 
previously determined. A sample of raw material with 
an exact weight of 300.0 grams was loaded into the 
mesh cassette of the extractor, after sealing the 


extractor, the unit was turned on. To determine the 
rate of loss of the residual oil content of the raw 
material in the first experiment, the plant was stopped 
after 100 seconds and after removal of the solvent (by 
dropping the pressure), the residual oil content of the 
sample was determined. In the second experiment, the 
installation with a new sample of raw materials of 
300.0 g stopped after 200 seconds, the residual oil 
content of the sample was determined, etc. up to 600 s 
extraction. Thus, the graphical dependence of the 
change in the content of extracted substances in the 
raw material on the extraction time, plotted in Fig. 11. 
The rate of mass transfer is high in the first 200 s of 
the process, since at that time oil is extracted on the 
surface layers of the raw material particles. 



—*— model ■ experiment 


Fig.ll. Variation in oil content of petals of grape seeds by extraction time 
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The first is a period of rapid extraction, during which 
the process proceeds at the fastest rate and the second 
is a slow extraction, the speed of which is much less. 
In the period of rapid extraction, which lasts for 600 s, 
up to 60% of the substances contained in the raw 
material are recovered. This period is characterized by 
the extraction of substances of more free molecular 
structures, and the particles of the raw material close 
to the surface. The second period, which lasts 3.000 
seconds, is characterized by the extraction of 
substances contained in the depth of microcapillaries 
and inside intact cells. During the second period, up to 
30% of the oil remaining in the raw material is 
recovered. 

The experimental experiments carried out confirm the 
adequacy of the obtained mathematical model of 
extraction of oil-bearing material in the apparatus 
with critical carbon dioxide. They will be used in 
future to find optimal solutions for the apparatus and 
technology of oil extraction using the critical state of 
carbon dioxide 
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